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Abstract 

This  paper  presents  a  lumped  element  model  of  a 
piezoelectric -driven  synthetic  jet  actuator.  A 
synthetic  jet,  also  known  as  a  zero  net  mass-flux 
device,  uses  a  vibrating  diaphragm  to  generate  an 
oscillatory  flow  through  a  small  orifice  or  slot.  In 
lumped  element  modeling  (LEM),  the  individual 
components  of  a  synthetic  jet  are  modeled  as 
elements  of  an  equivalent  electrical  circuit  using 
conjugate  power  variables.  The  frequency  response 
function  of  the  circuit  is  derived  to  obtain  an 
expression  for  Q„,„/VAC  ,  the  volume  flow  rate  per 
applied  voltage.  The  circuit  is  analyzed  to  provide 
physical  insight  into  the  dependence  of  the  device 
behavior  on  geometry  and  material  properties. 
Methods  to  estimate  the  model  parameters  are 
discussed,  and  experimental  verification  is  presented. 
In  addition,  the  model  is  used  to  estimate  the 
performance  of  two  prototypical  synthetic  jets,  and 
the  results  are  compared  with  experiment. 

1  Introduction 

Synthetic  jet  actuators  have  been  the  focus  of 
significant  research  activity  for  the  past  decade 
(Smith  and  Glezer  1998).  The  interest  in  synthetic 
jets  is  primarily  due  to  their  utility  in  flow  control 
applications,  such  as  separation  control,  mixing 
enhancement,  etc.  (Amitay  et  al.  1998;  Smith  et  al. 
1998;  Chen  et  al.  1999;  Honohan  et  al.  2000; 
Chatlynne  et  al.  2001). 

A  schematic  of  a  synthetic  jet  actuator  is  shown 
in  Figure  1.  A  typical  synthetic  jet,  also  known  as  a 
zero  net  mass-flux  device,  uses  a  vibrating  diaphragm 
to  drive  oscillatory  flow  through  a  small  orifice  or 
slot.  Although  there  is  no  source,  a  mean  jet  flow  is 
established  a  few  diameters  from  the  orifice  due  to 

£ 

Graduate  Student,  Student  Member  AIAA 

¥ 

Associate  Professor 

Associate  Professor,  Member  AIAA 

! 

Assistant  Professor,  Senior  Member  AIAA 

§ 

Assistant  Professor,  Member  AIAA 

Copyright  ©  2002  by  the  University  of  Florida.  Published  by 
the  American  Institute  of  Aeronautics  and  Astronautics,  Inc.  with 
permission. 


the  entrained  fluid  (in  the  absence  of  any  grazing 
external  flow). 

J  Net  Flow 


Orifice 


Oscillating  Piezo-Composite  Diaphragm 


Figure  1 :  Schematic  of  a  Synthetic  Jet. 

In  addition  to  studies  that  emphasize 
applications,  there  are  numerous  others  that  have 
concentrated  on  the  design,  visualization,  and/or 
measurements  of  synthetic  jets  (Crook  et  al.  1999; 
Chen  et  al.  2000;  Crook  and  Wood,  2001;  Gilarranz 
and  Rediniotis,  2001).  Furthermore,  several 
computational  studies  also  have  focused  on 
fundamental  aspects  of  these  devices  (Krai  et  al. 
1997;  Rizzeta  et  al.  1998;  Mallinson  et  al.  2000; 
Utturkar  et  al.  2002).  Crook  and  Wood  (2001) 
emphasize  the  importance  of  understanding  the 
scaling  and  operational  characteristics  of  a  synthetic 
jet.  Clearly,  this  information  is  required  for  a  user  to 
design  an  appropriate  device  for  a  particular 
application.  In  addition,  feedback  control 
applications  require  the  actuator  transfer  function  that 
relates  the  input  voltage  to  the  output  property  of 
interest  (e.g.,  volumetric  flow  rate)  in  the  control 
system. 

The  design  itself  represents  an  electro¬ 
mechanical-acoustic  coupled  system  with  frequency 
dependent  properties  determined  by  device 
dimensions  and  material  properties.  The  analysis  and 
design  of  coupled-domain  transducer  systems  are 
commonly  performed  using  lumped  element  models 
(Fisher  1955;  Hunt  1982;  Rossi  1988). 

The  main  assumption  employed  in  LEM  is  that 
the  characteristic  length  scales  of  the  governing 
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physical  phenomena  are  much  larger  than  the  largest 
geometric  dimension.  For  example,  in  an  acoustic 
system,  the  acoustic  wavelength  must  be  significantly 
larger  than  the  device  itself.  If  this  assumption  holds, 
then  the  governing  partial  differential  equations  for 
the  distributed  system  can  be  “lumped”  into  a  set  of 
coupled  ordinary  differential  equations.  This 
approach  provides  a  simple  method  to  estimate  the 
dynamic  response  of  a  synthetic  jet  for  design  and 
control-system  implementation. 

The  purpose  of  this  paper  is  to  rigorously  study 
the  application  of  lumped-element  modeling  to 
piezoelectric -driven  synthetic  jet  actuators.  To  the 
authors  knowledge,  this  represents  the  first 
application  of  LEM  to  piezoelectric -driven  synthetic 
jets.  McCormick  (2000)  employed  LEM  to  a 
speaker-driven  synthetic  jet,  while  Rathnasingham 
and  Breuer  (1997)  were  the  first  to  develop  a  low- 
order  model  of  a  synthetic  jet.  In  this  paper,  the 
various  lumped  elements  for  each  component  of  a 
synthetic  jet  are  theoretically  developed.  The 
resulting  equivalent  circuit  is  then  analyzed  to 
understand  the  effects  of  geometry  and  material 
properties  on  important  design  parameters,  such  as 
resonance  frequency  and  volume  displacement  per 
applied  voltage.  The  model  assumptions  and 
limitations  are  discussed,  along  with  the  results  of  an 
ongoing  experimental  study  designed  to  quantify  the 
validity  of  this  modeling  approach. 

2  Lumped  Element  Model 

In  LEM,  the  coupling  between  the  various 
energy  domains  is  realized  via  simplified  device 
representations  connected  to  form  equivalent  circuits. 
The  equivalent  circuit  model  is  constructed  by 
lumping  the  distributed  energy  storage  and 
dissipation  into  ideal  generalized  circuit  elements.  In 
this  electroacoustic  analogy,  pressure  and  voltage  are 
effort  variables,  while  current  and  volumetric  flow 
rate  are  flow  variables.  For  this  paper,  we  employ  an 
impedance  analogy,  in  which  elements  that  share  a 
common  effort  are  connected  in  parallel,  while  those 
sharing  a  common  flow  are  connected  in  series.  For 
a  synthetic  jet,  three  different  energy  domains  are 
involved:  electrical,  mechanical,  and  fluidic/acoustic. 

The  electromechanical  actuator  consists  of  a 
clamped  axisymmetric  PZT  patch  bonded  to  a  metal 
diaphragm.  The  composite  diaphragm  is  driven  into 
motion  via  an  applied  AC  voltage.  The  primary 
purpose  of  the  piezoelectric  diaphragm  is  to  produce 
large  volume  displacements  in  order  to  draw  fluid 
into  and  out  of  the  cavity.  This  represents  a 
conversion  from  the  mechanical  to  acousto-fluidic 
domain.  Consequently,  we  focus  on  the  first 
axisymmetric  vibration  mode  of  the  composite 
diaphragm  (see  Prasad  et  al.  2002  for  details  of  this 


analysis).  Specifically,  linear  composite  plate  theory 
is  used  to  obtain  the  short-circuit  pressure -deflection 
characteristics.  Then,  the  diaphragm  is  lumped  into 
an  equivalent  acoustic  mass  and  acoustic  compliance. 
The  former  represents  stored  kinetic  energy,  and  the 
latter  stored  potential  energy. 

In  general,  the  cavity  contains  a  compressible 
gas  that  stores  potential  energy  and  is  therefore 
modeled  as  an  acoustic  compliance.  Finally,  viscous 
effects  in  the  orifice  dissipate  a  portion  of  the  kinetic 
energy  stored  in  the  motion  of  the  oscillating  fluid 
mass.  Therefore,  there  will  be  an  effective  acoustic 
mass  and  acoustic  resistance  associated  with  the 
orifice  neck.  Flow  through  the  orifice  also  produces 
losses  associated  with  the  discharge  of  flow  from  the 
jet  exit.  In  addition,  for  low  operational  frequencies 
where  kao  =C0ao/c0<SCl  (where  ao  is  the  orifice 
radius  and  Co  is  the  speed  of  sound),  an  acoustic 
radiation  mass  also  be  added  if  the  orifice  is  ejecting 
into  a  semi -infinite  medium  (Rossi  1988). 


Figure  2:  Equivalent  circuit  representation  of  a 
piezoelectric-driven  synthetic  jet. 

The  equivalent  circuit  representation  for  the 
synthetic  jet  is  shown  in  Figure  2.  In  the  notation 
below,  the  first  subscript  denotes  the  domain  (e.g., 
“a”  for  acoustic  and  “e”  for  electric),  and  the  second 
subscript  describes  the  element  (e.g.,  “D”  for 
diaphragm).  In  the  electrical  domain,  C,  is  the 
blocked  electrical  capacitance  of  the  piezoelectric 
diaphragm  driven  by  an  AC  voltage  VAC  .  The  term 
blocked  is  used  since  it  is  the  impedance  seen  by  the 
source  when  the  diaphragm  motion  is  prevented. 

In  the  acoustic  domain,  C  „  and  M  are  the 
acoustic  compliance  and  mass  of  the  piezoceramic 
composite  diaphragm,  respectively.  M  may 

include  a  radiation  mass,  if  appropriate.  Although 
not  shown,  an  acoustic  resistance  could  also  be 
included  in  series  with  C  n  and  M  n  to  model 
structural  damping  effects.  CaC  is  the  acoustic 
compliance  of  the  cavity,  while  RM  and  MaN  are 
the  acoustic  resistance  and  mass  of  the  fluid  in  the 
neck,  respectively.  Finally,  R  is  the  resistance 

associated  with  the  orifice  discharge,  and  MaRad  is 
the  acoustic  radiation  mass  of  the  orifice. 
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In  this  paper,  we  assume  that  the  synthetic  jet 
exhausts  into  a  semi -infinite  ambient  air  medium,  and 
that  the  diaphragm  is  not  subject  to  a  mean 
differential  pressure.  If  necessary,  a  vent  channel  can 
be  used  to  equilibrate  the  mean  static  pressure  across 
the  diaphragm,  in  a  manner  similar  to  a  microphone 
(Sheplak  et  al.  1998).  For  simplicity,  we  will  assume 
that  there  is  no  grazing  flow  and  neglect  nonlinear 
amplitude  and  compressibility  effects  in  the  orifice. 

The  structure  of  the  equivalent  circuit  is 
explained  as  follows.  An  AC  voltage  V AC  is  applied 
across  the  piezoceramic  to  create  an  effective 
acoustic  pressure  that  drives  the  diaphragm  into 
motion.  This  represents  a  conversion  from  the 
electrical  to  the  acoustic  domain  and  is  accounted  for 
via  a  transformer  with  a  turns  ratio  §a  with  units  of 
[Pa/V].  An  ideal  transformer  (i.e.,  power  conserving) 
converts  energy  from  one  domain  to  another  and 
obeys  the  relations 

i  =  t,Q  andVAC=-p  {1} 


In  addition,  a  transformer  converts  an  electrical 
impedance  Ze  to  an  acoustic  impedance  Za  via 


P/  P/ 

_  ^ AC  _  /tyg  _  / Q  _  Za 

i  ~tQ~t2~V 


{2} 


The  motion  of  the  diaphragm  can  either 
compress  the  fluid  in  the  cavity  or  can  eject/ingest 
fluid  through  the  orifice.  Physically,  this  is 
represented  as  a  volume  velocity  divider, 
Q  =  Q  +  Q  .  The  goal  of  the  design  is  to  maximize 
the  magnitude  of  the  volume  flow  rate  through  the 
orifice  per  applied  voltage  | Qou,  jV AC  \ . 


3  Equivalent  Circuit  Model  Analysis 

Before  estimating  the  lumped  parameters 
defined  above,  it  is  instructive  to  analyze  the 
equivalent  circuit  to  obtain  the  frequency  response 
function  Qout  (s)/VAC  (s) ,  where  s  =  j(i> .  Using  Eq. 

{2},  the  transformer  can  be  eliminated  by  converting 
each  of  the  acoustic  impedances  to  their  electrical 
equivalent.  The  result  is  depicted  in  Figure  3,  where 


Zen(s)  =  - 


aD  T 


sC 


ZAs)  =  - 


sC,„ 


,  and 


{3} 


[C Kn_  +  Ra0  )+  S{MaN  +  MaRad  )] 


Figure  3:  Alternative  equivalent  circuit  model. 

Substituting  in  the  expressions  for 

ZeD,ZeC,  andZe0  and  grouping  powers  of  s  in  the 

numerator  and  denominator  results  in 

t,CflDs  {4} 

VAC(s)  a4s4  +a3s3  +a2s2  +a,i+l 


where 


~  (Ko  +  RaN)(CaD  +C<,c)’ 

a2  =  (M aRad  +  MaN)  {CaD  +CaC  )  +  MaDCaD  , 
a3  =  CaCMaDCaD  (  +  RaN  )  >  and 

«4  =CaCMaDCaD(MaR„d+MaN)- 

Although  this  expression  is  complex,  it  reveals 
some  important  features  without  having  to  estimate 
any  of  the  parameters  in  Eq.  {5}.  For  the  purposes  of 
this  discussion,  these  parameters  can  be  thought  of  as 
constants,  although  in  reality  some  are  likely  to 
exhibit  some  frequency  and  amplitude  dependence 
(i.e.,  due  to  nonlinear  effects).  For  a  DC  voltage 
( s  =  0 ),  the  volume  velocity  is  zero.  At  low 
frequencies  ( s  — >  0 ),  the  volume  velocity  is 
proportional  to  dVA(M  since  the  transduction  factor 
is  defined  as  d)  =  d  /C  _  ,  where  d  is  an  effective 
acoustic  piezoelectric  constant  defined  in  Eq.  {19}. 
This  result  emphasizes  the  need  to  optimize  the 
design  of  the  piezoceramic  composite  diaphragm 
(Prasad  et  al.  2002). 

At  high  frequencies  (  s  — >  °°  ),  we  find  that 


V AC  CacCaDMaD  {M aN  +  M aRad  )  S 

The  output  therefore  decreases  at  a  rate  of  60 
dB/decade  and  is  inversely  proportional  to  the  masses 
and  compliances  in  the  system. 

The  denominator  in  Eq.  {4}  is  a  4th-order 
polynomial  in  s,  indicating  two  resonance 
frequencies.  It  is  difficult  to  obtain  a  compact 
analytical  expression  for  the  two  resonance 
frequencies  (Fischer  1955).  However,  the  resonance 
frequencies,  fl2,  are  controlled  by  the  piezoelectric 
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diaphragm  natural  frequency  /  and  the  Helmholtz 
frequency  of  the  cavity  fc 


f d  ~ 


2n  vM  aDCaD 


f  =  — 

J  C  ^ 


!7} 


2n  V  ( MaN  +  MaRad )  Cal 


If  MaD  »  MaN  +  MaRad ,  then  the  natural  frequency  of 

the  synthetic  jet  actuator  equals  that  of  the 
diaphragm.  At  resonance,  the  response  is  limited  by 
the  flow  resistances  in  the  orifice  and  the  acoustic 
compliance  of  the  diaphragm 


Qm 


V4  r  C  n  R  n  +  R  M 

AC  aD  aO  aN 


{13} 


with  the  constraint  that 


=  {8} 

To  gain  physical  insight  into  the  behavior  of  the 
device,  three  important  cases  are  examined  below. 

C 

3.1  Case  I:  Incompressible  Limit  — ^ — >  0 

^aD 

Assuming  that  the  fluid  is  an  ideal  gas,  then  the 
acoustic  cavity  compliance,  which  is  analogous  to  the 
inverse  of  the  gas  spring  constant,  is  obtained  from 
the  cavity  volume  V0 ,  gas  density  p0,  and  the  speed 
of  sound  c0  via 


cflC  =  - 


Poco 


{9} 


In  practice,  CcICD—. >0  is  achieved  by  minimizing 
the  cavity  volume  or  operating  in  a  liquid  medium. 

Since  the  coefficients  a}  and  a4  in  Eq.  {5}  are 
both  proportional  to  CaC  ,  the  synthetic  jet  transfer 
function  reduces  to  the  2nd-order  system 


Qou,  CO 

^ AC  (S) 


a 


2 


2 


S 


{10} 


where  the  prime  denotes  the  limit  with 

cjcaD^  0. 

Eq.  {10}  is  written  in  the  form  of  a  canonical 
2nd-order  system  s 2  +  2^(0„i  +  co2 .  By  inspection,  we 

find  that  the  natural  frequency  and  damping  ratio  for 
the  incompressible  case  are  given  by 


co 

incomp 


1 

\CaDMaD 

i  i  M  aN  +M  aRad 

Mr, 

{11} 


and 


r  =1(R  +R  ) 

^  incomp  ^  ' 


C 


MaN+MaRad+MaD 


.{12} 


3.2 


Case  II:  Rigid  Diaphragm  Limit  — 1 —  — >  0 

^ aC 


As  described  in  Prasad  et  al.  (2002),  the  size  of 
the  piezoceramic  patch  is  not  negligible  compared  to 
the  metal  diaphragm  for  high  actuation  performance. 
In  this  case,  the  piezoceramic  composite  diaphragm 
cannot  be  accurately  modeled  as  a  homogeneous 
circular  plate.  Nonetheless,  assuming  that  the 
diaphragm  is  clamped,  the  acoustic  compliance  of  a 
homogeneous  clamped  circular  plate  provides  insight 
into  the  scaling  behavior  of  the  diaphragm 


nab  (l-v2) 
1 6  Eh' 


{14} 


where  a  is  the  radius,  E  is  the  elastic  modulus,  v  is 
Poisson’s  ratio,  and  h  is  the  thickness.  From  Eq. 
{14},  CaD  decreases  with  decreasing  thickness  ratio 
a/h  and  increasing  elastic  modulus. 

As  in  the  previous  case,  the  coefficients  a3  and 
a4  in  Eq.  {5}  are  zero,  and  the  synthetic  jet  transfer 
function  reduces  to  a  2d-order  system.  The  limit 
Cao/Cac  ~ ^0  leads  to  the  following  expressions  for 
the  natural  frequency,  damping  ratio,  and  response  at 
resonance: 


co 


J _ 

M aRad  )  CaC 


{15} 


^sm=-(Ro  +  RN) 


c, 


M  aRad  +MaD 


,  {16} 


and 


Qou,  1  da 


V,  r  C  r  R  n  +  R  „ 

AC  aC  aO  aN 


{17} 


In  this  case,  the  natural  frequency  of  the  jet 
corresponds  to  the  Helmholtz  frequency.  At 
resonance,  the  response  is  limited  by  the  flow 
resistances  in  the  orifice  and  the  acoustic  compliance 
of  the  cavity.  By  comparing  with  Eq.  {13},  the 
resonant  response  differs  for  these  cases  by  the  ratio 
of  the  acoustic  compliances. 
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3.3  Case  III:  Equal  Resonant  Frequencies/;^ 

One  other  case  of  interest  occurs  when  the  two 
natural  frequencies  are  equal  to  each  other,  /  =  f2 . 
However,  this  does  not  imply  that  fD  =  f  ,  as  seen 

in  Eq.  {8}.  It  can  be  shown  that  in  this  case,  the 
following  is  true: 


r  C 

aC  | 

c  „  c. 


-  +  2=  ■ 


M  „ 


M  ..  +M  „  . 

aN  aRad 


{18} 


The  cavity  volume  can  thus  be  adjusted  accordingly 
to  provide  a  single  dominant  peak  with  large  resonant 
response  |  Qom  ft AC  |  • 

4  Model  Parameter  Estimation 

In  this  section,  we  outline  the  methods  and 
assumptions  used  to  estimate  each  of  the  quantities  in 
Eq.  {5}.  The  interested  reader  is  referred  to  Rossi 
(1988)  and  Beranek  (1993)  for  details  of  the 
methodology. 

The  piezoelectric  diaphragm  vibrates  in  response 
to  both  an  applied  AC  voltage  and  oscillatory 
differential  pressure  according  to  the  relation 

Q  =  j®{daVAC  +  CaDP) .  {19} 


Here  du  =  \(Q/jaVAC  )|p=„  is  the  effective  acoustic 

piezoelectric  constant  that  relates  the  volume  velocity 
of  the  diaphragm  to  the  applied  voltage  VAC . 

CaD  ={Q/ jd)P)\VAc=o=(AVolume/P)\v^=o  relates  an 

applied  differential  pressure  to  the  volume  change  of 
the  diaphragm. 

The  vertical  deflection  w(r)  due  to  an  applied 
differential  pressure  is  also  lumped  into  an  equivalent 
acoustic  mass  M aD  by  equating  the  lumped  kinetic 
energy  of  the  vibrating  diaphragm  to  the  total  kinetic 
energy  using 

\^ftlw(ft2nrdr,  {20} 

2  o  2 


where  p'(r)  is  the  distributed  mass  per  unit  area,  Q 
is  the  net  volume  velocity  of  the  diaphragm,  and 
w(r)  =  jcovv(r)  is  the  distributed  vertical  velocity. 
All  of  these  parameters  are  calculated  via  linear 
composite  plate  theory  (see  Prasad  et  al.  2002  for 
details). 

The  acoustic  resistance  of  the  neck  is  obtained 
assuming  fully -developed  laminar  pipe  flow  in  the 
neck  of  length  L  and  radius  a0 
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where  (t  is  the  viscosity  of  the  fluid  and  Q  is  the 
volume  flow  rate  produced  by  the  pressure  A P  . 

Using  the  same  assumption  of  fully -developed 
pipe  flow,  the  acoustic  mass  in  the  neck  is  obtained 
from 
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As  mentioned  above,  the  acoustic  radiation  mass 
MaRa d  can  be  modeled,  to  first  order,  as  a  piston  in  an 
infinite  baffle  if  the  circular  orifice  is  mounted  in  a 
plate  that  is  much  larger  in  extent  than  the  orifice 
radius 
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The  acoustic  resistance  associated  with  the 
discharge  from  the  orifice  can  be  approximated  by 
modeling  the  orifice  as  a  generalized  Bernoulli  flow 
meter  (White  1979;  McCormick  2000): 


R„ 
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{25} 


where  u  is  the  mean  velocity,  and  KD  =0(1)  is  a 

nondimensional  loss  coefficient  that  is  a  function  of 
orifice  geometry,  Reynolds  number,  and  frequency. 
Note  that  R  „  is  a  function  of  the  volume  flow  rate 

aO 

Qoiil  through  the  orifice  and  thereby  represents  a 
nonlinear  resistance. 

As  mentioned  above,  a  transduction  factor  is 
required  to  move  from  the  electrical  to  the  acoustic 
domain.  This  factor  is  given  by 

4>fl=7T--  (26} 

'-'aD 

The  blocked  electrical  capacitance  Ceb  in  Figure 
2  is  related  to  the  free  electrical  capacitance  of  the 
piezoceramic  C0  =eAp  /  hp  by 
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v  C0CaD  J 

where  K 2  is  the  electroacoustic  coupling  factor,  e  is 
the  dielectric  constant,  Ap  is  the  piezoceramic  area, 
and  hp  is  the  thickness  of  the  piezoceramic. 

The  next  section  describes  experiments  designed 
to  isolate  and  measure  these  parameters.  This  allows 
us  to  evaluate  the  effectiveness  and  limitations  of  the 
lumped  element  model. 

5  Model  Validation 

A  modular  piezoelectric -driven  synthetic  jet  was 
constructed,  as  shown  in  Figure  4,  to  perform  a  series 
of  experiments  to  test  the  validity  tf  the  lumped 
element  model  parameters.  The  modular  design 
permits  a  parametric  study  of  the  cavity  volume, 
orifice  diameter  and  length,  and  piezoelectric 
diaphragm  diameter.  In  addition,  an  access  hole  is 
provided  for  a  microphone  to  monitor  the  fluctuating 
pressure  inside  the  cavity. 

Onbc*  pltu 


*■  •  j  .  * 

Top  plate  ioiy  pl*t«  Dtiptragm  Scttom 

Cap  plata  plate 

Aertit  kali 

lot 
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Figure  4:  Assembly  diagram  of  modular  synthetic 
jet. 


sinusoidal  excitation  voltage  at  /  =  100  Hz.  The 

measured  natural  frequency  of  the  diaphragm  was 
-2100  Hz  while  the  computed  resonance  frequency 
was  -2800  Hz.  The  non-zero  displacement  at 
r/a  ~  1  in  Figure  5(a)  reveals  a  compliant  boundary 
that  does  not  achieve  the  “clamped”  boundary 
condition  assumed  in  the  theory.  When  another 
diaphragm  was  tested  using  a  separate  mounting 
system  and  properly  clamped,  the  agreement  between 
theory  and  experiment  was  excellent,  as  shown  in 
Figure  5(b).  The  measured  and  predicted  natural 
frequencies  agreed  to  within  2%. 


Table  1:  Piezoceramic  Diaphragm  Details. 


Shim  (Brass) 

I 

n 

Elastic  Modulus  (Pa) 

8.963X1010 

Poisson’s  ratio 

0.324 

Density  (kg/m2) 

8700 

Thickness  (mm) 

0.20 

0.10 

Diameter  (mm) 

23.5 

37 

Piezoceramic  (PZT-5A) 

Elastic  Modulus  (Pa) 

6.1X1010 

Poisson’s  ratio 

0.33 

Density  (kg/m3) 

7500 

Thickness  (mm) 

0.11 

0.10 

Diameter  (mm) 

20.5 

25.0 

Dielectric  Permittivity  (F/m) 

8.85xl012 

d3|  (m/V) 

-2.7X1010 

CaD  (s2.m4/kg) 

6.53X1013 

2.23X10'1' 

MaD  (kg/m4) 

8.15xl03 

2.43X103 

<pa  (Pa/V) 

88.6 

15.1 

Future  testing  will  assess  the  severity  of 
nonlinear  effects  when  the  excitation  amplitude  is 
increased.  We  will  also  seek  to  directly  measure  the 
acoustic  compliance  of  the  diaphragm  CaD  in  a 
normal  incidence  plane-wave  tube. 


C*=C0(1-k2)=C0 


5.1  Piezoelectric  Transduction 

The  first  experiment  tested  the  linear  composite 
plate  theory  that  provides  estimates  for  C  D  ,  M  , 

and  da .  This  was  accomplished  by  measuring  the 
velocity  of  the  clamped  vibrating  diaphragm  (excited 
by  V AC )  using  a  scanning  laser  vibrometer  and 

integrating  the  velocity  to  displacement  in  the 
frequency  domain.  The  clamped  circular  diaphragm 
was  removed  from  the  synthetic  jet  apparatus  and 
mounted  on  an  optical  table.  The  test  was  also 
performed  in  a  vacuum  chamber  to  eliminate  fluid 
loading  effects. 

Table  1  provides  the  dimensions  and  material 
properties  of  the  diaphragms.  Figure  5(a)  shows  a 
comparison  between  the  predicted  and  measured 
response  of  the  piezoceramic  diaphragm  to  a 
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(b)  Clamped  boundary  case. 


Figure  5:  Comparison  between  predicted  and 
measured  response  of  piezoceramic  diaphragms  to  a 
sinusoidal  excitation  voltage  at  low  frecpiency. 
Theory  (  )  and  scanning  laser  vibrometer  in 

ambient  air  ( O,  )  and  in  a  vacuum  chamber  (  O). 

5.2  Cavity  Acoustic  Compliance 

The  value  for  the  cavity  acoustic  compliance 
CaC  is  obtained  from  Eq.  {9}.  The  cavity  volume 

can  be  calculated  from  the  geometry.  To  test  the 
theory,  the  orifice  was  replaced  with  a  solid  cap  to 
provide  a  closed  cavity  and  all  leaks  were  carefully 
minimized.  The  piezoceramic  was  then  driven  with  a 
sinusoid,  and  the  displacement  of  the  vibrating 
diaphragm  was  measured  with  a  laser  displacement 
sensor.  A  B&K  1/8  in.  microphone  then  measured 
the  fluctuating  pressure  in  the  closed  cavity. 


Figure  6:  Measured  acoustic  compliance  CaC  vs. 
frequency  in  closed  cavity  of  synthetic  jet. 


The  amplitude  of  the  sinusoid  was  adjusted  to 
avoid  harmonics  in  either  signal.  The  frequency 
response  function  between  the  pressure  and 
displacement  signal  at  the  diaphragm  center  w0  was 
used  to  measure  P/w0  at  several  frequencies  and 
calculate  CaC  ,  as  shown  in  Figure  6.  Using  the 
average  measured  value  of  9.89  MPa/m  and  the 
measured  mode  shape,  the  cavity  volume  was 


determined  to  be  2.6xl0'6  m3.  This  is  within  4%  of 
the  cavity  volume  calculated  from  the  geometric 
drawings. 

5.3  Acoustic  Mass  and  Resistance  in  Orifice 

The  flow  in  the  neck  of  the  orifice  is  modeled  as 
a  steady,  fully -developed  laminar  flow  in  a  circular 
duct  of  radius  a0  and  length  L .  The  resulting 
Poiseuille  flow  leads  to  the  expressions  for  acoustic 
resistance  and  acoustic  mass  defined  in  Eqs.  {21}  and 
{23},  respectively.  At  higher  frequencies,  the 
velocity  for  the  case  of  flow  in  a  circular  channel 
driven  by  an  oscillating  pressure  gradient  is  discussed 
in  White  (1974) 


(  A  ■ 
u(r,t)=  j- 


cop.,L 


1— 


.CO  r- 


-M 


.  co a: 


eim  ,  {28} 


where  /  is  a  Bessel  function  of  zero  order,  and  v  is 
the  kinematic  viscosity.  The  velocity  u  is 
proportional  to  the  pressure  gradient  and  inversely 
proportional  to  p()co  .  Furthermore,  the  velocity 
profile  is  characterized  by  the  Stokes  number 
St  =  ^co aljv  .  In  the  limit  of  St  0  ,  the  velocity 

profile  asymptotes  to  a  parabolic  shape  (i.e., 
Poiseuille  flow).  As  St  increases,  the  thickness  of 
the  Stokes  layers  decreases  below  a0 ,  leading  to  an 
inviscid  core  surrounded  by  a  viscous  annular  region. 


Figure  7:  Variation  in  velocity  profile  vs.  Stokes 
number  St  =  1,  10,  and  30  for  oscillatory  channel 
flow  in  a  circular  duct. 


These  velocity  profiles  can  be  integrated 
numerically  to  obtain  Qout  and  obtain  the  real 
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(resistive)  and  imaginary  (reactive)  parts  of  the 
acoustic  impedance  as  a  function  of  the  Stokes 
number 

^L  =  ZaN  =  RaN  +  jXaN  =  RaN  +  juMaN  .  {29} 

The  results  shown  in  Figure  8  reveal  that,  at  low 
frequencies,  the  acoustic  resistance  asymptotes  to  the 
steady  value  given  in  Eq.  {21}  and  increases 
gradually  with  frequency.  However,  the  acoustic 
mass  is  approximately  constant  with  frequency.  The 
data  in  Figure  8  can  be  used  to  provide  frequency 
dependent  estimates  for  the  acoustic  resistance  and 
mass  in  the  lumped  element  model. 

io'  i . . — — — — - — • — — —5  I 
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St;=«a{/v 

Figure  8:  Nondimensional  resistance  and  reactance 
vs.  St2  for  oscillatory  channel  flow  in  a  circular 
duct. 

It  should  also  be  noted  that,  depending  on  the 
aspect  ratio  L/a{)  of  the  orifice,  the  fully  developed 
assumption  may  not  be  valid.  Only  for  large  values 
of  L/a0  is  the  fully -developed  assumption  expected 
to  be  reasonable.  For  small  values  of  L/a0  ,  the 

orifice  loss  given  in  Eq.  {25}  is  expected  to 
dominate. 

The  models  discussed  in  this  section  are  simple 
and  neglect  potentially  significant  issues,  such  as 
nonlinear  effects  due  to  large  amplitude  pressure 
oscillations  in  the  cavity  (Ingard  1967)  and  transition 
to  turbulent  flow  and  compressibility  effects  in  the 
orifice.  The  grazing  flow  effects,  which  are  relevant 
when  the  synthetic  jet  interacts  with  a  boundary 
layer,  have  also  been  ignored  (Mittal  et  al.  2001). 

5.4  Comparison  between  Model  &  Experiment 

In  this  section,  the  lumped  element  model  is  used 
to  predict  the  performance  of  two  synthetic  jets  vs. 


frequency.  The  dimensions  of  the  device  are 
summarized  in  Tables  1  and  2. 

A  Dantec  two -component  Laser  Doppler 
Velocimeter  (LDV)  system  was  used  to  measure  the 
magnitude  of  the  peak  velocity  produced  by  the 
synthetic  jet  vs.  frequency.  The  synthetic  jet  was 
placed  inside  a  large  glass  fish  tank,  which  was  filled 
with  oil  seed  particles  and  closed.  The  probe  volume 
was  located  over  the  center  of  the  orifice  as  close  as 
possible  to  the  exit  plane.  The  system  was  operated 
in  backscatter  mode  and  resulted  in  a  probe-volume 
length  that  was  larger  than  the  diameter  of  the  orifice. 
Because  of  the  poor  spatial  resolution,  the  peak 
velocity  in  the  jet  was  measured  using  a  phase¬ 
averaging  scheme. 


Table  2: 

Synthetic  Jet  Details. 

Cavity: 

I 

n 

Volume  (m3) 

2.50x  IO"6 

5.00xl(T6 

Orifice: 

Radius  (mm) 

1.65 

0.84 

Length  (mm) 

1.65 

0.84 

60 
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Figure  9:  Comparison  between  the  lumped  element 
model  and  experiment  for  Case  I. 

Case  I  uses  a  smaller  piezoceramic  diaphragm 
and  cavity  volume  and  a  larger  orifice  than  Case  II. 
However,  both  devices  use  an  orifice  with  the  same 
aspect  ratio  L/a0  and  are  driven  with  VAC  =  50  V 
amplitude  sinusoids.  The  model  prediction  in  Figure 
9  has  the  correct  shape,  but  the  frequency  of  the  first 
resonance  peak  is  high.  Furthermore,  the  velocities 
are  overpredicted.  Note  that  the  peak  velocity  is 
estimated  from  Poiseuille  flow,  where 

Mmax  =  2H  =  2(20„,/(ft  a2 ) .  Since  this  relation  is  only 

valid  at  low  frequencies  for  ducts  with  fully - 
developed  flow,  this  approximation  is  likely  to 
introduce  significant  error  as  frequency  increases. 
Eq.  {28}  can  be  used  to  determine  the  relationship 
between  u  and  O  , ,  but  this  has  not  been  done  in 

max  zz'out  7 
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this  paper.  Clearly,  the  velocity  profile  must  be 
measured  and  integrated  to  obtain  Q  . 

Case  II  corresponds  approximately  to  the 
situation  described  in  Section  3.3.  In  particular,  this 
combination  brings  the  two  resonance  peaks  close 
together,  resulting  in  a  single  dominant  peak  with  jet 
velocities  of  over  60  m/s.  The  lumped  element 
model  accurately  predicts  the  resonance  frequency 
and  maximum  velocity  and  also  possesses  the  proper 
shape  of  the  frequency  response  function. 


Frequency  (Hz) 

Figure  10:  Comparison  between  the  lumped  element 
model  and  experiment  for  Case  II. 


These  results  reveal  the  power  and  shortcomings 
of  the  model  in  its  present  form.  Clearly,  the  flow  in 
the  vicinity  of  the  orifice  must  be  studied  further  in 
order  to  obtain  better  estimates  of  the  losses  and 
velocity  profile  characteristics.  Furthermore,  from 
Figure  8  we  expect  the  resistance  to  increase  with 
frequency,  but  the  current  model  only  uses  the  DC 
value  in  Eq.  {21}.  Furthermore,  the  loss  coefficient 
Kd  in  Eq.  {25}  is  treated  as  a  constant  but  is 
probably  a  function  of  (at  least)  the  orifice  Reynolds 
number,  Stokes  number,  and  geometry. 

6  Conclusions  and  Future  Work 

A  lumped  element  model  of  a  piezoelectric - 
driven  synthetic  jet  actuator  has  been  developed  and 
compared  with  experiment.  LEM  provides  a 
compact  analytical  model  and  valuable  physical 
insight  into  the  dependence  of  the  device  behavior  on 
geometry  and  material  properties.  The  model  reveals 
that  a  synthetic  jet  is  a  4th-order  coupled  oscillator. 
One  oscillator  is  a  Helmholtz  resonator,  and  the 
second  is  the  piezoelectric  diaphragm.  Simple 
arguments  reveal  three  important  special  cases 
corresponding  to  single  oscillators.  One  case  occurs 
when  the  fluid  is  incompressible  (i.e.  water)  or,  in  the 
case  of  a  gaseous  medium,  the  acoustic  compliance 
of  the  cavity  is  small  compared  to  that  of  the 
piezoelectric  diaphragm.  The  second  case  is  similar 


to  that  of  a  rigid  piston  and  occurs  when  the  acoustic 
compliance  of  the  piezoelectric  diaphragm  is  small 
compared  to  that  of  the  cavity.  In  this  case,  the 
synthetic  jet  acts  Ike  a  driven  Helmholtz  resonator. 
The  third  case  occurs  when  the  cavity  compliance  is 
adjusted  to  provide  a  single  resonant  peak  as 
described  in  Eq.  {18}. 

Methods  to  estimate  the  parameters  of  the 
lumped  element  model  were  discussed  in  some  detail 
and  experiments  were  performed  to  isolate  different 
components  of  the  model  and  evaluate  their 
suitability.  The  results  indicate  that  the  linear 
composite  plate  theory  is  accurate  when  the  assumed 
clamped  boundary  condition  is  achieved.  Similarly, 
the  cavity  acoustic  compliance  model  was  validated. 

The  details  of  the  flow  in  the  orifice  requires 
further  study.  It  is  this  region  that  dictates  the 
acoustic  mass  and  resistance  in  the  neck.  Accurate 
knowledge  of  the  acoustic  mass  is  required  to 
determine  the  Helmholtz  frequency  of  the  synthetic 
jet,  while  the  resistance  limits  the  response  of  the 
device  near  resonance. 

The  model  was  applied  to  two  prototypical 
synthetic  jets  and  found  to  provide  fair  agreement 
with  the  measured  performance.  Better  models  are 
needed  for  the  unsteady  flow  in  the  orifice,  including 
that  of  entrance  and  exit  effects.  In  addition, 
nonlinear  and  grazing  flow  effects  remain  to  be 
studied  in  a  rigorous  fashion. 

In  future  work,  additional  parameters  will  be 
varied  in  these  experiments  to  yield  a  performance 
database  for  the  synthetic  jet.  Emphasis  will  be 
placed  on  the  ratio  of  the  orifice  length  to  the  hole 
radius.  This  variable  will  be  systematically  varied  in 
concert  with  the  other  important  nondimensional 
parameters,  such  as  the  orifice  Reynolds  and  the 
Stokes  numbers.  Additional  velocity  measurements 
with  improved  spatial  resolution  will  also  be 
performed  to  map  out  the  spatial  variations  in  the 
synthetic  jet  velocity  field.  Finally,  an  optimization 
study  of  the  entire  lumped  element  model  is 
underway  and  should  lead  to  synthetic  jets  with 
improved  performance. 
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